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I. INTRODUCTION 

Entanglement is a key property that makes quantum 
information theory different from its classical counter- 
part 0]. Maximally entangled states, for example, are 
the basis of quantum state teleportation which has re- 
cently been demonstrated Under realistic condi- 
tions, however, one will only be able to generate par- 
tially entangled mixed states. It is then of interest to 
be able quantify the amount of entanglement in such 
states. Some measures of entanglement for mixed states 
have been suggested recently |^|-||. They are useful, for 
example, as upper bounds to entanglement purification 
protocols [p|-p^ and to the quantum channel capacity of 
certain quantum communication channels [pljslJl^. 

Unfortunately entanglement measures for mixed states 
(which are relevant in the presence of noise) are usually 
quite hard to calculate analytically, although an analytic 
expression for the entanglement of formation of two spin- 
1/2 particles is now known [|5| . The general case remains 
unsolved. For some problems, however, it is not so im- 
portant to know the exact amount of entanglement (a 
quantity that is not unique anyway). It would be com- 
pletely sufficient if one knew which state of a family of 
states has the most entanglement. To answer this ques- 
tion it would be sufficient to find a (hopefully as simple 
as possible) quantity that preserves the ordering of den- 
sity operators with respect to entanglement, i.e. that for 
two measures Ei and E2 and any two density operators 
pi and p2 we have that Ei{pi) > {<)Ei{p2) is equivalent 

to E2{pi) > {<)E2{p2)- 

In this paper we will compare the entanglement of for- 
mation for two spin-1/2 particles, for which a closed an- 
alytical form is known, with a quantity which was pro- 
posed in as a way to quantify the degree of entan- 
glement of a mixed state. The basis of this 'measure 
of entanglement' is the Peres-Horodecki criterion for the 
separability of bipartite systems [ p^jjl^ . Given a state 
of, for example, two spin-1/2 systems one calculates the 
partial transpose of the density operator. The state is 
separable exactly if the partial transpose is again a pos- 
itive operator. If, however, one of the eigenvalues of the 
partial transpose is negative then the state is entangled. 
One can now imagine that the amount of entanglement 



is quantified by the modulus of this negative eigenvalue, 
i.e. the larger it is, the larger the entanglement of the 
state. 

It is important to check whether this measure indeed 
preserves the ordering of density operators as it has been 
used for this purpose in some publications. 

In section II we will explain the entanglement of for- 
mation and subsequently we summarize some properties 
of the negative eigenvalue measure of entanglement. Fi- 
nally, we present in section III both a numerical and an 
analytical comparison of the two measures of entangle- 
ment with respect to the ordering of density operators 
induced by them. In section IV we sum up the results of 
this paper. 



II. ENTANGLEMENT MEASURES 

In this section we briefly describe some entanglement 
measures in particular the entanglement of formation and 
the negative eigenvalue measure of entanglement. 

There are not very many good measures of entangle- 
ment. One example is the relative entropy of entangle- 
ment [|-|. In fact, it gives rise to the most restrictive 
upper bound on the channel capacity of the depolariz- 
ing channel |||,|l^. Unfortunately, even for two spin-1/2 
particles no general analytical expression has been found 
for it so far although many special cases can be solved 
analytically. 

A second measure of entanglement and actually the 
first measure of entanglement that has been proposed for 
mixed states is the entanglement of formation . It ba- 
sically describes the amount of entanglement that needs 
to be shared previously in order to be able to create a 
particular ensemble in a given state p by local opera- 
tions. Mathematically, this means that we find that pure 
state ensemble that realizes the state p and which has 
the smallest amount of entanglement, i.e., 

Ef{p) = mm y2piE.vN{\^i){i}i\) , (1) 

where is a set of not necessarily orthogonal pure 

states. The entanglement of formation is known to be 
larger than the relative entropy of entanglement which 
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proves that in general quantum state purification meth- 
ods cannot recover all the entanglement that has been 
invested in the creation of the quantum state ||^. 

A nice feature of the entanglement of formation is the 
fact that it can be solved analytically for a system of 
two spin- 1/2 particles This allows for fast numeri- 
cal studies as the cumbersome minimization Eq. (|l|) can 
be avoided. 

The entanglement of formation can be expressed in 
terms of the function 



£{C) 



(2) 



where 



= — xlogj X — (1 — x) log2(l — 2^) . (3) 
For a density operator p one defines the spin flipped state 

P= [Py® Py)p*iPy(g) Py) (4) 

where the * denotes the complex conjugate in the stan- 
dard computational basis |00), |01), |10), |11). One then 
finds the entanglement of formation to be 

Epip) ^ £iC{p)) (5) 

where the so called concurrence is defined by 



C{p) = max{0, Ai - A2 - A3 - A4}. 



(6) 



Here, Ai,...,A4 are the eigenvalues, in decreasing order, 
of the Hermitean matrix R — \J ^[pp^fp- For properties 
of this measure of entanglement the reader should con- 
sult the literature It is interesting to note that 
since £ as a function of the concurrence C is a strictly 
monotonous function and maps the interval [0, 1] on [0, 1] 
C can in fact also be regarded as a measure for entangle- 
ment. 

We will now consider the negative eigenvalue of the 
partial transpose of a density operator as a measure of 
entanglement. In the next section we will then compare 
it to the entanglement of formation. 

For two spin-1/2 particles (which form the two systems 
A and B) any disentangled state p can be written as the 
convex sum of product states 



P - 



E 



(7) 



States which permit a representation of the form Eq. (j^) 
are also called separable. For two spin-1/2 particles there 
is a simple criterion to decide whether a given state is sep- 
arable or not 16|. One calculates the partial transpose 
P"^^ of the density operator p in the computational ba- 
sis. This means that we transpose only one subsystem, 
either subsystem A ox B. If the resulting matrix is posi- 
tive semidefinite then the density operator p is separable; 



otherwise it is not. Therefore the partial transpose of an 
entangled state has a negative eigenvalue and the idea of 
the negative eigenvalue measure is to use the modulus of 
the negative eigenvalue to quantify the entanglement of 
the state p. In a mathematical form this reads as 



En [p) = I min{0, Xl^ , A^- , AJ- , AJ- } | 



(8) 



where the Af^ are the eigenvalues of the partial transpose 
P"^^ . However, we do not know whether this way of quan- 
tifying the entanglement constitutes a proper measure of 
entanglement. Therefore we will investigate in the next 
section numerically whether the entanglement of forma- 
tion and the negative eigenvalue measure are compatible 
from a different point of view. 



III. ORDERING INDUCED BY 
ENTANGLEMENT MEASURES 

We would expect that any two 'good' entanglement 
measures should generate the same ordering of the den- 
sity operators. This means that for two entanglement 
measures Ei and E2 and any two density operators pi 
and p2 we have that 



Ei{pi) > Ei{p2) ^ E2{pi) > E2{P2) 



(9) 



Why do we expect this relation to be true? If in one 
measure of entanglement pi contains more entanglement 
than p2 then we would expect that a quantum state pu- 
rification method would generate more singlets from an 
ensemble in state pi than ensemble p2. If £'2 is also a 
measure of entanglement then we would expect that p2 
would yield more singlets than pi. While this reasoning 
is not strict, it nevertheless indicates that Eq. (^) should 
be true for two 'good' measures of entanglement. 

In the next two subsections we will now check whether 
Eq. is satisfied for the entanglement of formation Ep 
and the negative eigenvalue measure of entanglement En- 



A. An analytical comparison 

For some classes of density operators we can easily 
check analytically whether the relation Eq. (^) is true. 
For pure states it is sufficient to consider states of the 
form = ajOO) 4- /3|11) as it follows from the Schmidt 
decomposition The entanglement of formation Ep 
then reduces to the von Neumann entropy of entangle- 
ment while the negative eigenvalue measure yields 



(10) 



Both measures decrease monotonically with a so that for 
pure states Eq. (p|) is satisfied. It should be mentioned 
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that for pure states |V')('0l with = a|00) + the 
concurrence C is given by 



C{\il)){iP\) ^ 2ayJ\ - a2 



(11) 



hence, also for arbitrary pure states p — |V')(V'I the neg- 
ative eigenvalue measure and the concurrence are con- 
nected by the simple equation 



C{p) - 2E^{p). 
Werner states are defined as 



(12) 



(13) 



where = (|01) - |10))/V2 is the singlet state and 
F S [1/4, 1]. The concurrence of Werner states is given 
by C{pf) = — 1 for > 1/2 (otherwise is a sepa- 
rable state) while the entanglement of formation is found 
to be 

£^f(pf) =-Mlog2Ai-(l-M)log2(l-A*) (14) 



with /i = 1/2 + ^yF(Y^^F). The negative eigenvalue 
measure simply gives 



(15) 



which means that Eq. ( |12D is also valid for Werner states. 
Again, both measures decrease monotonically with F and 
therefore Eq. (^ is satisfied. 

Moreover, no counterexample to Eq. (^) can be con- 
structed from pure states and Werner states, i.e. no two 
density operators pi and p2 can be found which violate 
Eq. (H) when pi is assumed to be a pure state and p2 
corresponds to a Werner state. Because of the monotony 
of £ in Eq. (|) this can already be seen when comparing 
the concurrences of pi and p2 with the negative eigen- 
value measure of the respective states: from Eq. ( |l^ it 
then immediately follows that 



C{pi) > C{p2) A En{pi) < En{p2) 



(16) 



cannot be fulfilled for any value of F. 

One may therefore suspect that Eq. (H) also holds for 
arbitrary mixed states. To tackle the question whether 
Eq. (j|) is generally satisfied and hence the two entan- 
glement measures induce the same ordering we have em- 
ployed Monte Carlo simulations. If Eq. (^) is violated at 
all, it is - roughly speaking - also of interest to investigate 
how 'badly' it is violated after all. In the following sub- 
section we present the results of the numerical test that 
we have performed. 



B. A numerical comparison 

In order to test the ordering induced by the two entan- 
glement measures under consideration we have numeri- 
cally generated a million pairs of random density matri- 
ces with respect to a certain well defined distribution. In 
the case of two spin-1/2 particles, which is the case of in- 
terest in this paper, one has to generate random density 
operators - i.e. linear, self adjoint, positive semidcfinite 
operators p of finite trace {— 1) - acting on the four di- 
mensional Hilbert space isomorphic to . One can 
represent p by the decomposition 



4 

E 

1=1 



(17) 



with pairwise orthogonal projections, where pi > and 
'^iPi — 1- This in turn corresponds to a 4 x 4 matrix of 
the form 



P 



= UDU\ 



(18) 



where the diagonal matrix D is defined by Dij = SijPi. 

A plausible choice for the ensemble of unitary 4x4 ma- 
trices which can be used to construct density matrices ac- 
cording to Eq. (^8|) is the one with the normalized Haar 
measure on the group of unitary matrices U (4) which 
is called the circular unitary ensemble |T^. It has been 
shown in [ pO[ that random unitary matrices representa- 
tive for the circular unitary ensemble can be obtained as 
follows: Let us set 

U^U^''^\^12,A2,X12) (19) 

X C/('^') (^23, fe, 0);7(^''' (013, V'13, Xia) 

X U^'^'^Hhl, V'34, 0)U<^^'^^ (024, V'24, O)C/(l'4)(0i4, ^14, X14), 

where the complex 4x4 matrices U'^^'^\ i,j = 1,...,4, 
with three real parameters 4>, x, ^^'^ V' are given by 



1 

sm (pe 
cos (fie 
cos (pe 
— sin 0e 
0, 



k — i, I = j, 
k — I ~ i, 
k^l ^ j, 
k = j, l = i, 
otherwise. 



(20) 



We now take ipij and Xij to be independent random vari- 
ables with a uniform distribution in the interval [0, 27r), 
and for a given random number ^ distributed uniformly 
in the interval [0,1) we choose (/>ij to be arcsin(^^/'^^*)) 
for i ~ 1,...,3. Matrices generated in this way are 
then random unitary matrices of the wanted type (ex- 
cept for a random phase which is of no significance for 
our purposes). Finally, a set of random density matri- 
ces can be obtained by appropriately choosing the di- 
agonal matrices D. More precisely, the random vector 
with entries pi, ...,p4 should be uniformly distributed on 
the manifold defined by J^iPi = 1- According to |2|] 
this can be achieved, e.g., by setting pi — 1 — £,1^^, 
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P2 (1 - C2^^)(l - Pi), P3 = (1 - - Pi - P2), and 
P4 = 1 - Pi - P2 - P3, where again, ^1, ^2, and ^3 are 
random numbers drawn with respect to a uniform distri- 
bution in the interval [0, 1). 

To investigate whether Eq. is satisfied we now 
draw random density matrices from the previously de- 
scribed ensemble and calculate the eigenvalues of the par- 
tial transposes and the eigenvalues of the products of the 
respective density operators with the spin flipped states. 
After discarding those density matrices which correspond 
to separable states we then determine values for Eff and 
Ep as described in section II. Finally, we check for pairs 
(pii P2) whether the sign of 



and 



AEf{pi,P2) 



A£^Ar(pi,P2) 



Ef{pi) - Ef{p2) 
Ef{pi) + Ef{p2) 



EnIpi) - En{p2) 
EnIpi) + En{p2) 



(21) 



(22) 



is identical. Fig. |l] shows a diagram in which Ai?j?(pi, P2) 
is plotted versus AEff{pi, P2) for 10"* pairs (pi, P2) of ran- 
dom density matrices. Although there is obviously a cer- 
tain correlation between AEp and Ai?jv, there are dots 
in the second and the fourth quadrant of the diagram 
which are associated with pairs of states which do not 
satisfy Eq. (^). One can therefore conclude that in the 
case of two spin- 1/2 systems this relation is not satisfied 
for arbitrary mixed states, and hence, the ordering in- 
duced by the entanglement of formation Ep and the one 
induced by the negative eigenvalue measure E^ is not 
the same. 

AEf{pi,P2) 



AEn{pi,P2) 

FIG. 1. AEf{pi, P2) versus AEn{pi, P2) for 10'' pairs 
of entangled states {pi,p2). Pairs of states with 
AEf(pi, P2)AEm{pi, P2) < are represented by dots in the 
second and fourth quadrant. 



It is interesting to note that most of the randomly 
drawn density operators do not correspond to entangled 
but to separable states; so more density matrices have to 
be discarded than can be kept for the investigation. In 
fact, from the numerical simulations we can estimate the 
probability Pp that a created mixed state is entangled as 
Pe ~ 0.365 ±0.001, which is in complete agreement with 
the analytical and numerical findings presented in ]2l[ | . 

From the Monte Carlo simulation we can also give a 
rather accurate estimate of how likely it is to find a pair 
of states which violates Eq. (j^, given that both states 
are entangled. From the relative frequency of a violation 
for a million pairs of density matrices we obtain for the 
estimate of this probability Py 



Pv w 0.047 ±0.001. 



(23) 




FIG. 2. Distribution of states: En versus Ep- 

Now that we know that the ordering is different it is 
interesting to see to what extent the values of En{p) and 
Ep{p) differ for a given state p, or how these random 
states are distributed in a diagram where Ep is plotted 
versus E^. Fig. ^ shows such a plot. Since the order- 
ing induced by Ep and Ejsi is not the same, there is no 
strictly monotonous function / : [0,1] [0,1/2] such 
that En{p) — f{Ep{p)) for all states p. This is also 
obvious from Fig. ||. In Fig. ^ again the distribution 
of states is shown, but this time the negative eigenvalue 
measure is plotted versus the concurrence C. We can see 
that most of the dots are located close to the diagonal 
connecting (0,0) and (1/2,1); this diagonal corresponds 
to states satisfying Eq. (12). Note that the numerical 
simulation also strongly suggests that for a state p with 
a certain value of C{p) the upper bound for the possible 
values of the negative eigenvalue measure En (p) is given 
by C(p)/2, that is, that in general C{p) > 2Ei\f{p) holds 
for any state p. 
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Furthermore, since for pure states the ordering induced 
by the two entanglement measures is the same but in gen- 
eral it is not, it is of interest to investigate the dependence 
of the probability that Eq. (^) is violated for a pair of 
entangled states {pi, P2) on the 'purity of those states' 
with respect to a certain characterization of the purity. 
In Fig. ^ we show the relative number of pairs (pi, P2) of 
states which do not satisfy Eq. (||) versus S — Si + S2 
based on a million pairs of entangled states. Here, Si 
and 5*2 are the linear entropies 



Si = tr{p, - pf), i = 1,2, 



(24) 



of pi and p2, respectively, as commonly employed, e.g., 
in decoherence studies Obviously, pure states cor- 

respond to a vanishing linear entropy. 



E 



N 



c 

FIG. 3. Same as Fig. ^ with Epi versus C. 

We observe that the fraction of pairs {pi.p-i) of states 
with AEn{pi, P2)^Ef{pi, P2) < increases monotoni- 
cally with the sum of the linear entropies of the respec- 
tive states. This indicates that the more mixed the two 
states are with respect to the arithmetic mean of their 
linear entropies, the larger is the probability that this 
pair violates Eq. (^. Above a certain cut off no pairs of 
entangles states can be found at all - the numerical data 
shown in Fig. ^ give an estimate for the probability den- 
sity of finding a pair of entangled states with a certain 
value of S (compare also [^). 

It should finally be mentioned that in yet an- 
other measure of entanglement incorporating the partial 
transpose is proposed: For a given state p the quantity 
E = J2i=i \^J^ I — 1 is taken as a measure, where again, 
the Xj'^ denote the eigenvalues of the partial transpose 
p-^^ of p. While the actual value of E for a given state 
is of course different from that of Epf, from numerical 
simulations we come to the same conclusion that the or- 
dering induced by this measure is different from the one 
induced by the entanglement of formation Ep. 




FIG. 4. Estimate for the probability density of finding a 
pair of states {p\,p2) with a certain value of S = 5i -I- 5*2 
(+). The solid line is a cubic spline interpolation. Appar- 
ently, above a certain cut off there are no pairs of entangled 
states any more. The histogram shows the relative number 
of pairs violating Eq. (^ compared to the total number of 
pairs of entangled states for a given value of S. Note that the 
scaling is different for both plots: the right axis belongs to 
the estimate for the probability density (-I-), the left one to 
the relative number of violating pairs (histogram) . 



IV. SUMMARY 

We have compared the entanglement of formation with 
a potential measure of entanglement that is given by the 
negative eigenvalue of the partial transpose of the den- 
sity operator of the system. In particular the ordering of 
density operators with respect to the amount of entan- 
glement induced by the two measures has been compared 
both numerically and analytically. We have shown that 
the negative eigenvalue measure does not induce the same 
ordering as the entanglement of formation. Therefore we 
do not expect it to be a 'good' measure of entanglement. 
In particular it cannot, in general, be used to determine 
the most entangled state for a given family of density 
operators as it has been used previously. 
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